Stimulated Raman adiabatic passage in a A-system in the presence of quantum noise 
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We exploit a microscopically derived master equation for the study of STIRAP in the presence of 
decay from the auxiliary level toward the initial and final state, and compare our results with the 
predictions obtained from a phenomenological model previously used [P. A. Ivanov, N. V. Vitanov, 
and K. Bergmann, Phys. Rev. A 72, 053412 (2005)]. It is shown that our approach predicts a much 
higher efficiency. The effects of temperature are also taken into account, proving that in b-STIRAP 
thermal pumping can increase the efficiency of the population transfer. 
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I. INTRODUCTION 

Adiabatic theorem [l[ provides a very powerful tool for 
quantum state manipulation, and indeed it has been the 
basis of many applications, aimed at generation of quan- 
tum states 0] or at the realization of quantum gates Q 
based on geometric phase Even if the wide range of 
validity of the adiabatic theorem has been recently crit- 
icized very recently a sufficient condition that guar- 
antees adiabatic evolutions has been proven Q. 

Based on the adiabatic theorem, stimulated Raman 
adiabatic passage (STIRAP) i] allows the transfer 
of population from a quantum state of a physical sys- 
tem toward another state, through an auxiliary inter- 
mediate state Q. The passage occurs via a dark state, 
which aligns with the initial state in the beginning of the 
process, and then gradually changes its structure toward 
alignment with the target state. The process ends when 
the dark and target states coincide. For the opposite de- 
lay of the two couplings between the three states, and for 
nonzero single-photon detuning, it is also possible to re- 
alize the b-STIRAP process, which instead exploits the 
adiabatic change of a bright eigenstate of the Hamilto- 
nian from the initial state to the target state. The main 
difference between STIRAP and b-STIRAP is that in the 
latter case the auxiliary state is effectively involved in the 
dynamics, in the sense that a certain amount of popula- 
tion is temporarily transferred to it during the process. 
This circumstance makes the b-STIRAP more sensitive 
than STIRAP to the presence of decay from the auxil- 
iary level. In fact, while in the absence of environmental 
interaction and classical noise the transfer from the ini- 
tial state to the target state is predicted to be perfect, 
in the presence of dissipative dynamics the efficiency of 
the process is negatively affected. Many different mod- 
els have been considered to study the effect of dephasing 
[loj and spontaneous emission from the auxiliary level ei- 



ther toward external states or toward internal states 
[lH, i.e. the initial and the target state. In all these 
models, the incoherent dynamics has been taken into ac- 
count phenomenologically. Very recently, a microscopic 
model to describe the STIRAP and b-STIRAP processes 
under dissipation from the auxiliary state toward exter- 



nal states has been presented [l3| . The derivation of the 
master equation from a model of interaction between the 
three-state system and a bosonic environment, and the 
relevant dynamics, show some interesting deviations from 
the predictions related to the phenomenological counter- 
part. In particular, the microscopic model predicts a 
much higher efficiency in the STIRAP scheme. 

In this paper, by using the same rigorous microscopic 
model as described above [l]|, we investigate the effect 
of spontaneous decay from the intermediate state inside 
the A-system. Because the initial and final states in STI- 
RAP are usually ground or metastable, the intermediate 
state is necessarily an excited state, which may decay 
both inside and outside the system |bl| - tl8| . The loss of 
efficiency caused by external decay is more detrimental 
for it leads to irreversible population loss from the sys- 
tem; it is also easier to describe and understand [Tll[l3j. 
The effect of internal decay on STIRAP is a much more 
subtle effect because the loss of efficiency is compensated 
by the concomitant optical pumping 12] . Here we de- 
velop a rigorous microscopic theory of internal decay in 
STIRAP and b-STIRAP, which reveals some unexpected 
features compared to the phenomenological model [l2| ■ 

Starting from a microscopic model of system- 
environment interaction, we derive a time-dependent 
master equation that describes the dynamics of our sys- 
tem. Then, after considering the resolution of the master 
equation, we compare the predictions coming from our 
model with the results coming from the phenomenolog- 
ical description of an analogous decay scheme [l2|. We 
find that the efficiency of the STIRAP process is higher 
than predicted before. Exploiting the master equation 
at nonzero temperature, we also study the effects of tem- 
perature, showing that the thermal pumping dramati- 
cally and negatively affects the efficiency of the popula- 
tion transfer in the STIRAP process, while has a slightly 
positive effect in b-STIRAP. 

The paper is organized as follows. In the next section 
we present the derivation of the Markovian master equa- 
tion of a system with time- independent Hamiltonian and 
a time-dependent system-environment interaction term. 
In the third section we apply the result of the previ- 
ous section and the general theory of Davies and Spohn 
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[13 . |20( to derive the master equation of out three-state 
system. Then, in the fourth section we show the results 
obtained at zero temperature and compare them with 
the results coming from the phenomenological model. In 
section V we show the effects of temperature, and finally 
in the last section we give some conclusive remarks. 

II. GENERAL FORMALISM 

In this section we consider the general problem of 
the derivation of the master equation for a system 
whose Hamiltonian H s is constant, while the system- 
environment interaction Hamiltonian contains oscillation 
terms. This will be useful in the next section, where 
we deal with a system described (in a rotating frame) 
by a slowly-varying Hamiltonian and interacting with a 
thermal bath through oscillating terms. According to 
the general theory by Davies and Spohn [2(|, un- 
der the hypothesis that the environmental correlation 
time is much smaller than the timescale of the Hamilto- 
nian change, we can treat this system by assuming that 
the system Hamiltonian is time-independent during the 
derivation of the master equation, and putting the time- 
dependence of the jump operators only after the deriva- 
tion. 

Therefore we start by considering a time-independent 
system Hamiltonian H s and the following time- 
dependent system-bath interaction Hamiltonian: 

H i = Yl ( A « elUa 1 + A « e ~ iUa ') ® B * ■ w 

Following the approach presented in [22[, let us intro- 
duce, for each Bohr frequency u> 

At(u)= J2 m^tW), (2) 

e' — e— U 

where n(e) is the projector on the subspace of the system 
Hilbert space corresponding to the energy eigenvalue e 



and the sum is extended over all the couples of energies 
e and e' such that e' — e = u>. The operators defined in 
this way satisfy both 

[H at At(u,)] = -u a A±(w), (3) 

and 

(A±{u,))* = AZ(-w), (4) 

giving 

e iH '*A±(u)e- iH - t = e- iut A±(u>), (5) 

e lH ° l (A±( w )) t e- lH ° l = e iut (A±( w )) t . (6) 

Another important property is that summing over all 
the Bohr frequencies (both positive and negative) one 
rcobtains the initial operators: 

At=J2 A H")- CO 

UJ 

In the Schrodinger picture we thus have: 
Hi = * + A « H e ~ tUla l )®B a , (8) 

which in the interaction picture with respect to H s + Hb 
becomes: 

Hj = e _i " * {A+ («) ' + A~ (u) *) ® B a (t) , 

a.cu 

(9) 

or, taking the Hermitian conjugate: 

Hi = (KM)* e-^ + {A-{^ e *--*) 

® Bi(t). (10) 
The formal resolution of the Liouville equation gives: 



d.str s {Hj(t - s)p(t)p B H I (t) - Hj{t)Hj(t - s)p(t)p B } + h.c. 



(11) 



from which, substituting the expansions of Hj, one gets the following master equation: 
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EE^^'^r^H (A+{u,)p(A+(w')? - (A+(^A+(u)p 

^^ e i(.-.'-a,,- Wa )t r +- (w) ^- (w)/3 (^ (w '))t _ (^+( w '))t A - (w)p 



/l.C. 



(12) 



r 



with 



and 



dse ^-^)«(st(t)B^(t- s )) ,(13) 



dse^^^BKiJB^i-j)) .(14) 



This is the most general form of the Born-Markov 
master equation before a Rotating Wave Approxima- 
tion (RWA) is performed. Under the hypothesis that 
uj a ,ujp 3> u),cj', one can single out very clear conditions 
for RWA. The only terms which survive are those for 
which uj a and ujp appear in the combination tu a — ujp 
with a = j3, and u = ui': 

p = ee r ~H ( A Mp( A a(^y 

uj a 

(A+iurf A+(u) P ) 

uj a. 

- {AM) 1 A-{uj)p)+h.c. (15) 



which, neglecting the Lamb shifts and coming back to 
the Schrodinger picture, becomes: 



uj a 

l{(A+(u>)yA+(uj),p} 

uj a. 

l{{A-(w)?A-(w),p}\ , 



(16) 



where 7 ++H = 2SR{r++M} and vM 
25R{r-"H}. 

III. OUR MODEL 



A. The system 



We consider a three-level system in A-configuration 
whose Hamiltonian is: 



Hsys(t) - 



fip(t)e- l (" 21 




-A)t 



f2 p (t)e 

n,(t) 



z(aJ2i— A)t 
„i(w23-A)t 





-i(w 2 3-A)t 



(17) 



w 3 



The system interacts with a bosonic bath. The free while the system-bath interaction Hamiltonian is 
bath is described by 

H mt = E.9i 12) (l 1 )( 2 l + I 2 )( 1 l)(^ + & D 



H B =J2^kblh, (18) 



E<7f ) (|3>(2| + |2)<3|)(^ + &i). (19) 
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FIG. 1: (Color online). States |1) and |3) are coherently cou- 
pled to level 1 2} (red bold lines). The state |2) is coupled to 
the other two states by a dipolar system-environment inter- 
action (green dashed lines). 



In the rotating frame associated with the transforma- 
tion T{t) = e iuJlt |1) (1| + e^ 2 "^' |2) (2| + e iUst |3) (3|, 
the total Hamiltonian: 



H = H s (t) + H B + H T (t) , 



with 



H a (t) 



o n p (t) o 
n p (t) a n s (t) 
o n s (t) o 



(20) 



(21) 



and 



fl>(t) = Y,9k 2) (e l(wi -" 2+A)t |1) (2| 

k 

+ e -*(«™+A)t |2)(1|) (b k + b{) 

+ J29k 2) (e i( " 3 -" 2+A)t |3) (2| 
fe 

+ e -*(<«s-e*H-A)t |2)(3|) (bk+bl). (22) 

The eigenstates and eigenvalues of H s (t) are: 
cj + = Qq cot ip, = 0, cj_ = — f2gtan(p, (23a) 



|+) = sin sin |1) + cosip |2) + sin</?cos0 |3) , (23b) 
|0) =cos0|l) -sin0|3) , (23c) 
| — ) = cos <p sin 1 1) — sin |2) + cos tp cos |3) , (23d) 



where 



tan0(t) 
tan 2</?(t) 



n p (t) 
n.(t) 

2fi(t) 



A(t) ' 

n(t) = ^/fip(i) 2 + n s (t) s 



(24a) 
(24b) 
(24c) 



It is well known 0] that, for the intuitive sequence of 
pulses, in which the probe pulse fl p (t) precedes the Stokes 
pulse one has | — ) = |1) for t — > — oo and | — ) = |3) 

for t — >• oo. Therefore, if all the pulses vary adiabatically, 
the population of the state |1) can be transferred to the 
state |3): this process is called b-STIRAP. On the other 
hand, if f2 s (£) precedes p (i), one has |0) = |1) for t — > 
—oo and |0) = |3) for t — > oo: in this counterintuitive 
sequence the population from |1) to |3) is adiabatically 
transferred through the state |0) and the process is called 
STIRAP. 



B. Master Equation 



Since in the rotating frame we have a slowly-varying 
system Hamiltonian and a time-dependent system- 
environment interaction term, we can use the general 
formalism presented before in order to derive the mas- 
ter equation that describes the dynamics of our system. 
We obtain the following master equation (for details see 
appendix [A]) : 



5 



= -i[H a ,p]+ [7+ + (^+o)cos 2 0cosV + 7 h + 6 + (^ + o)sin 2 0cosV] (\0) (+1 P 1+) <0| - \ {|+> (+| , p}j 

+ [7«."(<Oo-) cos 2 0sin 2 93 + 7» - (wo_) sin 2 sin 2 <0| p |0) {-| - i {|0) (0| , p} j 

+ [ita ) sin2 cos 4 93 + 7~ (w + _) sin 2 (9 sin 4 99 + 7++ cos 2 cos 4 93 + 7bb ~ (w+_) cos 2 (9 sin 4 93] 

x f|-> <+|p|+> <-|-|-CI+> <+|,p}) 
+ [(7+ + (0) + 7„T (0)) sin 2 sin 2 93 cos 2 <p + ( 7 ++(0) + (0)) cos 2 sin 2 93 cos 2 <p] 

x (1+) (+1 - 1-) H)p(l+> (+1 - h) H) - ^ {l+> (+1 + 1-) H,p) 

+ [7 < T."(wo+) cos 2 cos 2 93 + 7»"(wo+) sin 2 cos 2 93] (|+) (0| p |0) (+| - \ {|0> (0| , p}) 

+ [ 7a + +(^-o) cos 2 flsin 2 93 + 7b + h + (^-o) sin 2 flsin 2 93] (jO) <-| p |-> (0| - 1 {|-> (-| , p}) 
+ [7aa"( w -+) sin 2 (9 cos 4 93 + 7 ++(w-+) sin 2 6> sin 4 93 + 7 ~(u;_ + ) cos 2 6> cos 4 93 + 7 ++(w_ + ) cos 2 (9 sin 4 93] 
x (\+){-\p\-){+\-\{\-){-\,P}) , 



(25) 



where Ulnm = — w TO . 

From (fT5|) and Q14p one gets that the decay rates are 
given by a spectral density Jj{u>) multiplied by a factor 
depending on the photon population N(oS) of the bath 
modes at the relevant frequency corrected with ±u; a or 
ioJb depending on the case, i.e.: 



where g a (ui) (g&(w)) is the system- reservoir coupling con- 
stant <?fc 12 ^ (fffe 32 ^) in the continuum limit, and d(ui) is the 
reservoir density of states at frequency ui. 



7++(o;) = J^u - u)j)(l + N(uj - Uj )) 

Iji M = + "iX 1 + N ( u} + ^o)) 
^ + Ju) = J j (\uj-cj j \)N(\uj-u j \) 

Ijj {w) = Ji(\u+Ui\)N(\u + u j \) 



with j = a,b and 



w 6 



oj\ - oj 2 + A , 
= cj 3 - w 2 + A . 



U)-U). 

id + id. 

UJ — LO^ 
UJ + UJ. 



>o, 

>0, 
<0, 
< 0, 



(27a) 
(27b) 



The zero temperature spectral density Jj{uj) for gen- 
eral bosonic reservoir is given by [2l|, [22( : 



Jj(u) = d(uj) \gj(u)\ 



(28) 



It is important to note that, under the hypothesis that 
'j 3> uj for any Bohr frequency uj between the dressed 
(2®Jates in the rotating frame (which is the usual case since 
uij are optical frequencies associated with the atomic 
transitions while id's are of the order of magnitude of 
the coupling terms Q's), and taking into account that 
the frequencies in ()27a[) are negative, the only condition 
satisfied are id — idj > and uj + uij < 0. Therefore, 
in (|2l))) , only the rates of the first and fourth classes are 
possible. Moreover, at zero temperature only the rates 
of the first class survive, since the number of photons in 
the reservoir is zero. In such a case the master equation 
becomes: 
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= -i[H 3 ,p}+ [7+ a + (^ + o)co S 2 0cosV + 75 + h + (^+o)sin 2 ^cosV] (|0> (+| p R (0| - \ {{+) (+| , p}^ 

+ [ 7a + +(w-o) cos 2 0sin 2 <p + J+ + (oj. q ) sin 2 6 sin 2 f\0) (-\p\-) <0| - 1 {|-> (-| , p}j 

+ [7+ + (^+-)sin 2 0cosV + 7 ++( W+ _)cos 2 0cosV] x (j-) <+| p |+> <-| - \ {|+> (+| , p}) 

+ [7a + a + (^- + )sni 2 ^mV + 7 h + b + (^- + )cos 2 0sinV] x (-| P R (+1 " 5 {h> H , />}) 

+ [ita (0) sin2 61 sin 2 <p cos 2 </? + 7^+ (0) cos 2 6 sin 2 v? cos 2 tp] 

1+) (+1 - h) (-I)MR <+| - I-) (-1) - J {R (+| + I-) (-1 , P } 



(29) 



+ 







FIG. 2: (Color online). Scheme of the decays for the dressed 
states. There are transitions from |+) to | — ) and vice versa. 
Both |+) and |— ) decay toward |0). The dephasing between 
the states +) and | — ) is not represented. 



This equation shows that at zero temperature there 
are the following processes: transitions from |+) to |— ) 
and vice versa, transitions from |+) to |0) and from — ) 
to |0), and a dephasing process involving levels |+) and 
— ) (see figure [5]). This suggests the idea that the damp- 
ing can help to transfer population to level |0), so that 
the efficiency of the counterintuitive sequence should be 
positively affected by the dissipation. 



IV. ANALYSIS OF THE EFFICIENCY AT ZERO 
TEMPERATURE 



We consider Gaussian laser pulses: 
fit 



^0 -(t-T/2f/T 2 
2 



-(t+r/2) 2 /T 2 



(30a) 
(30b) 



taking into account that we have the so called intuitive se- 
quence (which corresponds to b-STIRAP) when f2 p = O2 
and f2 s = fii, while we get the counterintuitive sequence 
(which corresponds to STIRAP) when il p = f)i and 

O s = f?2- 

The phenomenological model which we will compare 
with our microscopic model corresponds to the following 
master equation: 



-i[H,p]-^D, 



(31) 



with 



D = 



(Ti + r 3 )p2i 




(IV 
2(Ti 

(rv 



- 1^3)^22 

~T 3 )P32 





(Ti + r 3 )p 2 3 

-2r 3/ 9 2 2 



(32) 

which describes spontaneous emission from level 2 to lev- 
els 1 and 3 with rates Ti and r 3 , respectively. Such a 
master equation is related to the bare states and then 
turns out to be time-independent. 

Concerning the microscopic model, we assume fiat 
spectrum for both the transitions 2 — t 1, correspond- 
ing to J a (w), and 2—5-3, corresponding to J&(w). In 
particular we will assume J a (uj) = T and Jt(w) = aT. 
This may come, for instance, from the assumption that 
the dipole moments between the states |1) and |2) and 
between the states |3) and |2) are proportional, so that 
g^ 12 ' = ag^ 2 ^ for any k in cq. (fT9| . 



In this section we analyze the efficiency of both STI- 
RAP and b-STIRAP processes, by numerically studying 
the post-pulse population of the target state |3) and com- 
pare the prediction of our model with the predictions of 
a phenomenological model introduced in ref. [l2j ■ 



A. The counterintuitive sequence 

We first analyze the counterintuitive sequence, where 
the population is carried by the dark state |0). 
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P 3 (co) 




io" 2 10" 1 1 10 io 2 io 3 io 4 

r 

FIG. 3: (Color online). Counterintuitive sequence. Final 
population vs T (in units of T _1 and in logarithmic scale) 
according to microscopic (dashed blue line) and phenomeno- 
logical (solid red line) model. The relevant parameters are 
flo = 25 T _1 , r = 1.5 T~\ T A = 1, a = 1. 



P 3 (oo) i 




0.5 1 1.5 2 



r 

FIG. 5: (Color online). Intuitive sequence. Final population 
vs r (in units of T -1 ) according to microscopic (dashed blue 
line) and phenomenological (solid red line) model. The rel- 
evant parameters are f2o = 25T" 1 , r = 1.5T -1 , TA = 1, 
a — 1. The two curves are essentially coincident. 




FIG. 4: (Color online). Counterintuitive sequence. Final 
population vs both F (in units of T _1 and in logarithmic 
scale) and a according to the microscopic model. The relevant 
parameters are f2 = 25T -1 , r = 1.5 T" 1 , T A = 1. 



Figure [3] shows the comparison between the micro- 
scopic and the phenomenological models, with a = 1 
and ri = T2 = r. It is evident that the microscopic 
model predicts a very high efficiency (essentially one) for 
a wider range of T. This can be explained on the basis 
of the decay scheme in fig. [2j all the decay processes 
describe either jumps toward |0) or toward states which 
in turn decay toward |0), so that the dark state is robust 
against zero temperature dissipation. 

The robustness of the counterintuitive scheme is not re- 
lated to the special choice a = 1. Indeed, figure |4] shows 
the dependence of the post-pulse population of state |3) 
on both r and the parameter a which characterizes the 
difference in the intensities of the dipolar coupling con- 
stant involving different couples of levels. It is quite ev- 
ident that the efficiency of the scheme is not affected by 
a discrepancy in the decay rates. 



B. The intuitive sequence 

Concerning the b-STIRAP process (i.e. in the intu- 
itive sequence), we find that the two models predict very 
similar results. In particular, fig. [5] shows the perfect co- 
incidence of the predictions of the two models in the case 
a = 1 and ri = T2 = T. Moreover, from this figure one 
can see that (for both models) the efficiency is very sen- 
sitive to the presence of decays, so that it almost drops 
to zero at FT = 1 . The reason of the fragility of the ef- 
ficiency in this scheme is that, while all the populations 
are guided by the decay towards the state |0), population 
transfer is instead carried on by the state — ). 

It is worth noting that the result of the comparison is 
qualitatively quite similar to the result of the comparison 
we made in connection to the scheme with external de- 
cay [Hj]. Indeed, in both cases the predictions from the 
microscopic and the phenomenological models are almost 
coincident for the intuitive sequence, while for the coun- 
terintuitive sequence we find that the microscopic model 
predicts a higher efficiency. Nevertheless, we stress here 
the fact that the enhancement of efficiency in STIRAP in 
the strong damping limit is due to very different mech- 
anisms in the two cases of external and internal decay. 
In fact, while for external decay a very strong damp- 
ing is responsible for a dynamical decoupling of the dark 
state, which then is protected against losses, in the case 
of internal decay the dissipation is instead responsible for 
transitions toward the state that carries the population, 
therefore protecting the process of population transfer. 

V. ANALYSIS OF THE EFFICIENCY AT 
NONZERO TEMPERATURE 

In this section we consider the effects of nonzero tem- 
perature. Looking at (|26|). we see that the N(ujys are 
evaluated at very close frequencies, which are essentially 
lo « UJ21 ~ W23. For this reason, we described tempera- 
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ture by a single number TV which is the number of pho- 
tons in the reservoir modes of frequencies close to the 
bare atomic transitions. 

We have seen that at zero temperature, for increasing 
r (and even for quite small values of the decay constant) 
the efficiency falls down to zero when the intuitive se- 
quence is applied. When temperature is non-vanishing 
all the transitions included in (|25p but not present in 
(|29|) must be considered. In particular, transitions from 
0) to |— ) should increase the efficiency of the b-STIRAP 
process, since in this scheme the population is transferred 
through the state |— }. On the other hand, in the coun- 
terintuitive sequence we should get a lower transfer ef- 
ficiency, since thermal terms are responsible for loss of 
population from state |0) during the STIRAP process. 

Figure^ shows the dependence of the post-pulse pop- 
ulation of the state |3) on T and temperature (through 
the number of photons N in the relevant reservoir 
modes), for the intuitive sequence. It is well visible that 
the efficiency, which goes to zero for large T in the zero 
temperature regime, reaches nonzero values for non van- 
ishing temperature. Such efficiency reaches a maximum 
value for intermediate values of temperature. As an ex- 
ample, figure shows the temperature dependence (in 
a wider range) of the efficiency for T = 1: in this case 
the optimal point is reached at N ~ 10. 

Figure [7] shows the post-pulse population of the state 
|3) on r and temperature for the counterintuitive se- 
quence. In this case it is well visible that the temper- 
ature negatively affects the efficiency of the population 
transfer. Indeed, even a very small amount of thermal 
photons is responsible for a significant diminishing of the 
post-pulse population, which instead approaches unity at 
rigorously zero temperature. 



VI. DISCUSSION AND CONCLUSIVE 
REMARKS 
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FIG. 6: (Color online). Intuitive sequence, (a) Final popula- 
tion vs r (in units of T~ ) and the number of photons N. (b) 
Final population vs the number of photons N (in logarithmic 
scale) , for Y — 1. In both cases, the relevant parameters are 
Qo = 25 T -1 , t = 1.5 T -1 , TA = 1, a = 1. 



In this paper we have analyzed the effects on a STI- 
RAP scheme of the losses of the auxiliary level towards 
the two metastablc states, examined by means of the nu- 
merical resolution of a master equation which has been 
microscopically derived. It is worth noting that, in order 
to remove the rapid oscillations in the system Hamilto- 
nian, we are forced to describe the system in a rotating 
frame, where the system-bath interaction term turns out 
to be time-dependent. Therefore, when deriving the mas- 
ter equation we have to deal both with a slowly varying 
system Hamiltonian, which can be treated following the 
general theory by Davies and Spohn, and a rapidly os- 
cillating system-environment interaction term. This last 
point makes the final master equation different from what 
one usually gets. Indeed the zero temperature transi- 
tions are not guiding the system towards the dressed (ro- 
tating) ground state |— ), but instead towards the state 
|0). In this sense, the oscillating terms in the system- 
environment interaction act like a pumping which makes 




FIG. 7: (Color online). Counterintuitive sequence. Final 
population vs T (in units of T" 1 and in logarithmic scale) 
and the number of photons N (in logarithmic scale). The 
relevant parameters are f2o = 25 T -1 , r = 1.5 T~ , T A = 1, 
a = 1. 
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the counterintuitive sequence much more robust than the 
intuitive sequence. 

The inclusion of the nonzero temperature terms in 
the master equation partially modifies these conclusions. 
In fact, thermal photons switch on different transitions, 
which make the post-pulse population of state |— ) dif- 
ferent from zero. As a consequence the efficiency of the 
counterintuitive sequence is reduced, while on the other 
hand the efficiency of the intuitive sequence increases. 
However, such an increasing is not high enough to make 
the intuitive sequence preferable to the counterintuitive 
one. 

Though the analysis of this STIRAP scheme (with the 
same decay channels), performed in previous papers by 
means of phenomenological dissipative terms, has shown 
a high robustness of the counterintuitive sequence with 
respect to losses, our analysis shows that the scheme at 
zero temperature is much more robust for very high decay 
rates. This is due to the fact that, as already pointed out, 
the zero temperature decay channels force the system 
to the state |0), which is the population carrier in this 
sequence. 
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Appendix A: Derivation of the Master Equation 
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